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The calculated specific heat of optimally K-doped BaFe 2 As 2 in density functional theory is about 
five times smaller than that found in the experiment. We report that by adjusting the potential on 
the iron atom to be slightly more repulsive for electrons improves the calculated heat capacity as 
well as the structural, magnetic, and electronic properties of Bao. 6 Ko. 4 Fe 2 As 2 . Applying the same 
correction to the antiferromagnetic state, we find that the electron-phonon coupling is strongly 
enhanced. 

PACS numbers: 71.15.Mb, 71.20.-b, 74.70.Xa 


I. INTRODUCTION 

The discovery^ of superconductivity in LaFeAsOi_ x F x 
with a transition temperature of 26 K in 2008 triggered 
unprecedented interest and further research in iron-based 
superconductors. So far, superconductivity was found 
in four main families of iron-based compounds: 1111, 
122, 111, and 11^ These iron-based materials have two 
phases in the normal state: one is a paramagnetic metal 
and the other is an antiferromagnetic metal. Supercon¬ 
ductivity emerges in both the paramagnetic and antifer¬ 
romagnetic metal phases via application of hydrostatic 
pressure or carrier doping of the parent materials. Hence, 
it is expected that understanding the electronic and mag¬ 
netic structures of the metallic normal states of these sys¬ 
tems is a needed ingredient for unraveling the origin of 
the superconductivity of iron-based materials. 

Many experimental and theoretical studies have been 
done on the normal states of iron-based superconductors, 
and a consensus^ has been reached in these systems 
that the Coulomb interaction among the electrons is not 
strong enough to induce a Mott insulating phase. How¬ 
ever, the Coulomb interaction plays an important role in 
determining the electric and magnetic properties. In the 
early stages of this research, theoretical insight into the 
properties of these materials was gained by calculations 
based on density functional theory (DFT) within the lo¬ 
cal density approximation (LDA) or generalized gradient 
approximation (GGA). However, LDA and GGA have 
some limitations in describing the normal states of iron- 
based superconductors. In a paramagnetic phase, the 
measured mass of low-energy quasiparticles is 2-3 times 
larger than that calculated within LDA or GGA. In ad¬ 
dition, the measured magnitude of the ordered moment 
in an antiferromagnetic phase is 2-3 times smaller than 
that obtained with LDA or GGA. Furthermore, LDA 
and GGA studies related to the specific heat of these 
materials are not consistent with the experimental data. 
The theoretical Sommerfeld coefficienlpU— of optimally K- 
doped BaFe 2 As 2 is about five times smaller than that 
found in the experimental data*2r— . 

There have been many attempts to describe electronic 
correlations in these materials by combining LDA or 


GGA calculations with a dynamical mean-field theory 
(DMFT), quasiparticle self-consistent GW (QSGW), or 
the Gutzwiller methodJA— Using these methods, many 
of the electric and magnetic properties of correlated iron- 
based superconductors can be reproduced. For example, 
effective masses and Fermi surfaces (FSs) across all fam¬ 
ilies of iron compounds are well described in the frame¬ 
work of DFT+DMFT— and QSGW-i^, as well as or¬ 
dered moments and the fluctuations of local moments 
within DFT+DMFTi 14 i 16 However, a calculation of the 
electron-phonon coupling coefficient (which is needed for 
heat capacity and superconductivity estimates) within 
these approaches is non-trivial. Therefore, we use a sim¬ 
pler method to calculate electronic and magnetic prop¬ 
erties of these materials. 

In this work we study the heat capacity of 
Bao. 6 Ko. 4 Fe 2 As 2 superconductor (T c =38 within 

a semi-empirically modified GGA potential, following 
studies^ of an FeSe monolayer. We show that one can 
choose a small repulsive potential located on the iron 
atoms (+A term) so that the calculated specific heat coef¬ 
ficient is increased from 7„=12 mJ mol -1 K -2 in GGA to 
j n =38 mJ mol -1 K -2 in GGA+A, much closer to recent 
experimental findings (7^=40-50 mJ mol -1 K~ 2 )i2r— 
The increase in j n relative to GGA comes mostly from 
the increased density of states (DOS) at the Fermi level 
and to a smaller extent from an enhanced electron- 
phonon coupling. Since Bao. 6 Ko. 4 Fe 2 As 2 is near a mag¬ 
netic phase transition, we also computed the heat capac¬ 
ity in the striped antiferromagnetic ground state, present 
in the parent compound. Just as in the nonmagnetic 
calculation, we again find an increased j n (from 6 to 
12 mJ mol -1 K -2 ) when +A term is added. However, 
unlike in the nonmagnetic calculation, the increased y n 
originates mostly from increase in the electron-phonon 
coefficient A. 


II. METHODS 

Our calculations are based on ab-initio norm- 
conserving pseudopotentials and the Perdew-Burke- 
Ernzerho^S functional as implemented in the SIESTA 
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codei 21 Electronic wavefunctions are expanded with 
pseudoatomic orbitals (double-^ polarization). We treat 
the potassium doping within the virtual crystal approxi¬ 
mation. 

Following Ref. 0 we modify the GGA potential 
Vgga (r) by adding a repulsive potential on each iron 
atom in the calculation, 

Vgga(» + A^/dr-r*!). (1) 

i 

Here /(r) is a positive dimensionless function peaked at 
the nucleus of the iron atoms (r;) and the extent of /(r) is 
comparable with the size of d orbitals in the iron atoms. 
We discuss the choice of A and f(r) in subsection III Al 
The GGA+A approach can be understood as a variant 
of the constrained DFT (CDFT) formalism— The CDFT 
approach adds a general constraint to the density, 

! w c( r )p a ( r ) dr = N c ( 2 ) 

a d 

where w c ( r) acts as a weight function that defines the 
constrained property. The Kohn-Sham total energy is 
minimized under the constraint from Eq. [2] by making 
the following functional stationary, 


W\p, V c } = E\p\ + V c(Y,J <(rK(r)dr - N^J . (3) 


Here V c is a Lagrange multiplier corresponding to the 
constraint. Therefore, in the effective Hamiltonian of 
the CDFT formalism, there is an additional potential 
V c Wc (r) coming from the constraint. Since the GGA+A 
potential (Eq. [l]) has the same form as the constraint 
potential in the CDFT approach, GGA+A method has 
the same effect as constraining the number of electrons 
around the iron atom. 


A. Choice of Af(r) term 

Now we discuss our choice of the correction term Af(r ) 
appearing in Eq. [T] 

Following previous work on the FeSe monolayer— we 
first choose /(r) = e~ 1 ' 0r in atomic units (Bohr ra¬ 
dius) with the extend comparable with the size of the 
iron atom d-orbital. Second, we tune A from 0 up to 
A c until one of the properties of Bao. 6 Ko. 4 Fe 2 As 2 agrees 
better with the experimental data. We choose to tune 
the occupied bandwidth of the M-point electron pocket 
since it is severely overestimated in GGA (it is 130 meV 
in GGA while ~ 0-10 meV in the experimen t 23 ’ 24 ). We 
find that using A c = 1.3 Ry has the desired effect of 
tuning the M-point bandwidth to about 4 meV. 

Just as in Ref. |T^ we find that the choice of f(r) is 
not very important for physical properties as long as it 
is localized on the iron atom and A is tuned for each 


TABLE 1. A comparison of the arsenic height, iron-arsenic- 
iron angle, magnetic moment ( pp e ) on iron atom, and the 
energy difference (A E) per one iron atom between antifer¬ 
romagnetic stripe and nonmagnetic ground state in GGA, 
GGA+A, and from experiment (Ref. [jj]) in Bao.6Ko.4Fe2As2- 



As height 

(A) 

Fe-As-Fe 

n 

/i Fe 

(Mb) 

A E 
(eV) 

GGA+A 

1.44 

105.9 

2.26 

-0.19 

GGA 

1.30 

112.5 

2.91 

-0.33 

Experiment 

1.37 

109.9 




choice of /(r). For example, using Af(r) = 2.2e 17r 
or Af(r ) = 5.5e _3 ' 5r results in nearly indistinguishable 
band structure of Bao. 6 Ko. 4 Fe 2 As 2 . 

Using A = A c improves not only the occupied band¬ 
width of the M-point electron pocket but other proper¬ 
ties of Bao.eKo. 4 Fe 2 As 2 as well. For example, structural 
parameters relevant for superconductivity (arsenic height 
and iron-arsenic-iron angle )2£r— are both moved in the 
direction towards experimental value. Going from GGA 
to GGA+A the arsenic height is increased from 1.30 A 
to 1.44 A while the iron-arsenic-iron angle is decreased 
from 112.5° to 105.9°. In addition, antiferromagnetic 
ground state is suppressed in GGA+A. See Table [T] for 
more details. We confirmed that the tendency for im¬ 
proved structural and magnetic properties is independent 
of the choice of /(r). 

Finally, using A = A c the calculated heat capacity 
of Bao. 6 Ko. 4 Fe 2 As 2 is more than three times larger as 
compared to the GGA value, and in good agreement with 
the experimental value. We discuss heat capacity in more 
detail in Sec. hyi 


III. ELECTRONIC STRUCTURE 

Now we discuss the electronic structure of 
Bao. 6 Ko. 4 Fe 2 As 2 in GGA and GGA+A. In all of 
our calculations we perform a full structural relaxation 
for both forces on atoms and stresses on the cell. We 
sample the Brillouin zone on a uniform 32x32x32 
k-point mesh. 

Figure [T] compares the electron density in 
Bao. 6 Ko. 4 Fe 2 As 2 in GGA and GGA+A. From panel b of 
the figure it is clear that including the +A term transfers 
some of the electronic density from the iron atom to 
the outer region. The maximal change in the electronic 
density is about 7 % and it occurs on a charge density 
peak near the iron atom. 

Figure [5] compares the band structure and the Fermi 
surface in GGA and GGA+A. We compare these results 
to the experiment in Sec. IIII Al 

In the GGA case, as in a previous calculation,— there 
are three hole pockets at the zone center (T), and two 
electron pockets at the zone corner (M). However, the 
band structures and the Fermi surfaces in GGA+A are 
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FIG. 1. The electron density of Bao.6Ko.4Fe2As2 in the non¬ 
magnetic state on a line between the nearest-neighbor (left) 
and next-nearest-neighbor (right) iron atoms within GGA 
(blue) and GGA+A (red). Difference between GGA and 
GGA+A is shown in panel (b). Densities of both semi-core 
(3s, 3p) and valence (3d, 4s) states are included in our pseu¬ 
dopotential calculation. 


both quantitatively and qualitatively different in sev¬ 
eral respects. First, the occupied bandwidth of the M- 
point d yz and d zx electron pockets in GGA+A is 4 meV 
[Fig. [21(a)], while it is 130 meV in GGA. In addition, the 
effective mass of these pockets is increased by a factor 
of 3-4 in GGA+A and the shape of the Fermi pocket in 
GGA+A is more elongated towards the T and Z points. 

Second, the area of the hole pockets at T and Z is 
changed in GGA+A. Specifically, in GGA+A the size of 
the d xy hole pockets at T and Z is increased by a factor 
of 4, so that it is larger than remaining two pockets. 
In addition, the d z 2 hole pocket is not present at Z in 
GGA+A so that now there are only two hole pockets at 
the Z point (versus three hole pockets at Z in GGA). 
Therefore, a three-dimensional ellipsoidal Fermi surface 
exists at F in GGA+A. 
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FIG. 2. Electronic band structures and Fermi surface of 
Bao.6Ko.4Fe2As2 in the nonmagnetic states both in GGA+A 
(panels a and b) and in GGA (c,d). Dominant orbital char¬ 
acters (defined in the single-iron unit cell) are represented in 
blue {d xy ), red {d yz ), green {dzx), black {d z 2 ), and yellow 
(both d yz and d zx ) color. High symmetry points in the Bril- 
louin are defined in the two-iron unit cell. Reciprocal space 
axes fci and &2 are perpendicular to the tetragonal c-axis. 


served in angle-resolved photoemission spectroscopy from 
Ref. Third, the presence of three-dimensional FS in 
Bao. 6 Ko. 4 Fe 2 As 2 was suggested from c-axis polarized op¬ 
tical measurements^ The optical experiment found that 
the c-axis data only exhibit a small difference across T c . 
This indicates the existence of three-dimensional FS with 
a dispersive band along the c axis. 


IV. SPECIFIC HEAT 


A. Comparison with ARPES 


In this section we discuss the calculated specific heat 
of Bao. 6 Ko. 4 Fe 2 As 2 . The specific heat coefficient is 
defined as, 


Now we compare modifications in the band structure 
due to +A term with the currently available experimental 
data on Bao. 6 Ko. 4 Fe 2 As 2 band structure. 

First, in angle-resolved photoemission spectroscopy 
(ARPES) experiment, three hole pockets are observed at 
the zone center and the largest pocket is shown to orig¬ 
inate from d xy orbital^ as in our GGA+A calculation. 
Second, large elongation of M point pocket towards the 
r and Z points we find using +A was experimentally ob- 


7„ = (1 + A)7o. (4) 

Here 70 is the Sommerfeld coefficient proportional to 
DOS at the Fermi energy, and A is the electron-phonon 
coupling coefficient. 

First we discuss the density of states in GGA and 
GGA+A. In GGA the DOS at the Fermi energy of non¬ 
magnetic Bao. 6 Ko. 4 Fe 2 As 2 is 4.4 states eV -1 f.u. -1 (the 
energy dependence of DOS is shown in Fig. [3[). Similar 
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FIG. 3. Bao. 6 Ko. 4 Fe 2 As 2 density of states for both spin com¬ 
ponents, per two-iron formula unit in the nonmagnetic (NM), 
checkerboard antiferromagnetic, and single-stripe antiferro¬ 
magnetic state, both in GGA+A (a,b,c) and GGA (d,e,f). 


value (3.1-5.5 states eV 1 f.u. for DOS was found in 
previous calculations^ 

In our GGA+A calculation, DOS at the peak value 
near the Fermi level is 11.9 states eV -1 f.u. -1 , almost 
three times larger than in GGA. Since 70 is proportional 
to DOS, it is also increased by a factor of 3 in GGA+A 
over GGA (see Table HT1). 

Increase in the DOS after inclusion of +A term origi¬ 
nates from the changes of the band structure at the M 
point. In the GGA+A, the bottom of the electron-like 
band at the M point (and the corresponding van Hove 
singularity) is placed almost at the Ep. Furthermore, 
the DOS at the van Hove singularity is enhanced due to 
the renormalization of the band width and the formation 
of a saddle point at the M point [Fig. 0(a)]. 

After having discussed the 70 , we now discuss the con¬ 
tribution of the electron-phonon coupling coefficient (A) 
to the heat capacity . We calculated the electron- 
phonon coupling coefficient A using the Wannier in¬ 
terpolation technique^ and the Quantum-ESPRESSO 
packaged 3 The electron-phonon coupling in the nonmag¬ 
netic GGA+A calculation is 0.37, about two times larger 
than 0.18 obtained in GGA (see Table HU). However, the 
heat capacity ( 7 „) is proportional to 1 + A so the increase 
in A in GGA+A increases 7 „ by 16 %, in addition to the 
dominant increase from larger DOS. 

Taking both terms together (70 and 1 + A) we find 
that within GGA+A method specific heat coefficient 7 „ 
equals 38 mJ mol -1 K -2 , which is much closer to the 
experimentally measured values (40-50 mJ mol -1 K -2 ) 
than the GGA result (12 mJ mol -1 K -2 ). 


TABLE II. A comparison of the Sommerfeld coefficient ( 70 ), 
electron-phonon coupling (A), and enhanced normal-state 
specific heat coefficient (y„) of Bao. 6 Ko. 4 Fe 2 As 2 in GGA, 
GGA+A, and from experiment (Ref. Isl-flDl. Coefficients 70 
and 7 „ are given in mJ mol -1 K -2 in a two-iron atom unit 
cell and for both spin components. 



70 

A 

7n 

Nonmagnetic 

GGA+A 

28 

0.37 

38 

GGA 

Experiment 

10 

0.18^- 

12 

40-50 

Checkerboard 

GGA+A 

4.1 

0.80 

7.3 

GGA 

26 

0.3321 

34 

Single-stripe 

GGA+A 

6.6 

0.90 

12 

GGA 

5.1 

0.1822 

6.0 


A. Antiferromagnetic ground states 

So far we discussed the specific heat in the non- 
mangetic ground state of Bao. 6 Ko. 4 Fe 2 As 2 , now we con¬ 
sider two antiferromagnetic ground states: striped and 
checkerboard. The striped case is especially important, 
since this is the experimentally determined ground state 
of the parent compound BaFe 2 As 2 - We study the alter¬ 
native ground state (checkerboard) for a comparison with 
the striped phase. 

The (single-)stripe order consists of ferromagnetically 
arranged chains of iron atoms, with antiferromagnetic 
aligment between neighboring chains. On the other hand, 
in the checkerboard antiferromagnetic order magnetic 
moments on all neighboring iron atoms in point in op¬ 
posite directions. 

For easier comparison with the nonmagnetic calcula¬ 
tions, in our magnetic GGA+A calculations we use the 
same value of A c and the same function f(r). 

In the striped state, the peak in DOS occuring 50 meV 
below the Fermi level in GGA is shifted to Ep —230 meV 
when +A is included but there is no significant change 
in the DOS at the Ep [Figs. [3]+) and (f)]. However, 
in the checkerboard state within GGA+A we obtain the 
DOS at Ep equal to 1.7 states eV -1 f.u. -1 , which is 
about one sixth of the GGA result (see Figs. [3)(b) and 
(e)). This suppression in the checkerboard state is due 
to the occurrence of a Jahn-Teller distortion in GGA+A, 
which is lowering the crystal symmetry from tetragonal 
to orthorhombic. 

Eventhough within GGA+A DOS at Ep is relatively 
small in the striped state (2.8 states eV -1 f.u. -1 ) the 
electron-phonon coupling is significantly larger than in 
the nonmagnetic case. We obtained A = 0.90 (see Ta¬ 
ble ED in striped state which is ~60 % larger than in 
GGA. As we said earlier, DOS in striped state is nearly 
the same in GGA and GGA+A. Therefore, strong en¬ 
hancement of A in GGA+A must originate from other 
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sources, and not simply from increased DOS. However, 
the origin of this enhancement is not the focus of this 
paper, and it will be reported elsewhere. 

V. CONCLUSION 

By increasing the potential on iron atoms, making 
them slightly more repulsive for electrons significantly 
improves the structural, magnetic, and electronic prop¬ 
erties ofBa 0 . 6 Ko. 4 Fe 2 As 2 , as calculated within DFT. The 
main result of this paper is that with a corrected potential 
(+A) on iron atom, the heat capacity of Bao. 6 Ko. 4 Fe 2 As 2 
is increased more than threefold, in good agreement with 
experimental data. Applying the same correction to the 
magnetic states, we find that electron-phonon coupling 
is strongly enhanced. This observation might be crucial 


in understanding the superconducting properties of iron- 
based superconductors. 
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